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Abstract 

Recurrence relations derived via the Chetyrkin-Tkachov method 
of integration by parts are apphed to reduce scalar three-loop bub- 
ble (vacuum) diagrams with a mass to a limited number of master 
integrals. The reduction is implemented as a package of computer 
programs for analytic evaluation in FORM. The algorithms are appli- 
cable to diagrams with any integer powers on the lines in an arbitrary 
dimension. A physical application is the evaluation of the three-loop 
QCD correction to the electroweak rho parameter. 
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Vacuum (bubble) Feynman integrals (without external momenta) appear 
as low-energy limits of certain physical amplitudes or as Taylor coefficients 
of multipoint Green functions. The coefficients can then be used to recover 
the functions in the whole complex plane of momenta The presence of 
virtual heavy particles, like the top quark in the quantum chromo dynamics 
(QCD), generates an effective high-energy scale, which makes perturbation 
theory applicable owing to asymptotic freedom. In quadratically and linearly 
divergent diagrams the contributions of heavy particles are enhanced by the 
power of the mass, so that light particles can well be considered as massless. 
This leaves us with an important special case of just one nonzero mass. 

After performing the Dirac and Lorentz algebra, any vacuum diagram 
can be reduced to some linear combinations of scalar bubble integrals. In 
the three-loop case with the full tetrahedron topology of the diagram, all 
scalar products in the numerator can be expressed through the quadratic 
combinations in the denominators. Thus, in the most general case, only 
a product of some powers of the denominators should be integrated. A 
prototype defines the arrangement of massive and massless lines in a diagram. 
Individual integrals are specified by the powers of the denominators, called 
indices of the lines. 

In the dimensional regularization with N = 4 — 2e, any massless bubbles 
are trivially equal to zero, so that at least one massive line should be present. 
Fig.|l| displays all possible three- loop prototypes. It was convenient also to 
distinguish some reduced prototypes with a line missing, -E2-4, since they are 
generated in evaluating several different full prototypes. The Bm and Bj^ 
types have been completely analyzed in Ref. and need not be discussed 
here. 

The method of recurrence relations [§, 0] connects integrals of the same 
prototype but with different values of the indices. Using these relations 
ingeniously enough, one can reduce any integral to a limited set of so-called 
master integrals. However, that remains still a kind of art without any strict 
assertions as to the minimal set of the master integrals or the most efficient 
strategy. Let us derive a relation ||^] for a generic triangle subgraph with 
masses on its lines mi, 1712, and ms, line momenta pi, P2 = Pi — P12, and 
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Figure 1: The three- loop scalar bubble prototypes with one mass. Double 
(single) lines refer to massive (massless) propagators in the momentum repre- 
sentation. Any line may have an integer index, a power of the denominator. 
The indices stay as arguments of the corresponding functions. Numbers de- 
fine the ordering of the arguments. 
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P2. = Pi- Pi3, and indices ji, j2, and ja, respectively: 

d pi f d^pi / .... 2m^ 



J dpi 4' 4' 4' J clVsVg^ ^ 



Cl 



, + m2 - + Ci2 - Ci , m^ + m^- m^g + C13 - Ci n 
+ J2 +J3 ), (1) 

C2 C3 ^ 

where = pl + ml. Dividing or multiplying by just increments or decre- 
ments index jk of the k th line. The corresponding operator can be denoted 
by [0. For an arbitrary L-loop diagram, integrating by parts the first 
derivatives provides us with linearly independent relations in total. It is 
convenient to refer to triangle recurrence relations by specifying only the line 
numbers: {123} for Eq. (|1|), line 1 being the base. Half sum of the relations 
for three faces of a tetrahedron as their base lines form a triangle, like {124}, 
{534}, and {623} in D^, gives a relation {dim} that is evident on dimensional 
grounds: 

[| iV - ji - ... - j6 + m\j, 4+ + j5 5+ + J6 6+)]D3(ji, Je) = . (2) 

Now follows a brief description of evaluating various prototypes. In Fig. |l| 
they have been ordered from the 'most difficult' to the 'simplest'. Whenever 
in Dq a denominator is absent {jk < 0) the diagram is immediately reduced 
to the 'simpler' type by expanding the power of the polynomial. A typical 
trick to bring any positive index down to 1 is as follows. A combination is 
sought, like 3{146} — {416} — {614} for ji > 1 in Dq, in which only one 
'highest' term m^jil"*" is present, all others having the sum of the indices less 
by one. The highest term can be expressed through the others until the index 
on the line reaches 1. Thus we arrive at the master integral Dq{1, 

If D5 has js < 0, it is reduced to Bj^; ji_4 < to D4. As jg > 0, it is 
profitable to solve {612} (no for a massless exchange between particles of 
unchanged masses Q]) with respect to the free term. Eventually, this brings 
J3,4 or Je to zero. To get rid of the numerator jg < 0, a denominator on an 
adjacent line is typically used: if for example ji > 1, {315} can be solved 
relative to l+G". Otherwise, as ji_4 = 1, the quadratic denominators can 
be created by solving {126} with respect to the free term. The irreducible 
master integral is .05(1, 1, 1, 1, 1, 0). 

In D4, je < leads to Bm] < to Dm; iz,A < to D^] and ji,2 = 
to £'4. The numerator on linel can be eliminated by {246} if jg > Ij by 
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2{dim}-{534} if jg > 1; by {dim}-{345} if > 1; by {dim}-{624} if 
j2 7^ 1; or by {215} otherwise. Solving {215} relative to m^l+ diminishes 
the denominator ji > 1; {512} reduces js > 1; {435} — {534} reduces > 1; 
{136}+{246} reduces jg > 1, leading to 0^(1, 1). 

As j4 6 < in D3, this is D2; j5 < ^ Dj^; j2 < ^ Bj^. The numerator 
ji < is normally reduced by 2{dim} — {534} as > 1; by {623} as j2 > 1; 
by {236} as Jq > 1; by 2{dim} — {435} as J5 > 1; and {dim} creates J^^q > 1 
otherwise. The special case ji = js = is more efficiently worked out as E4 
with js = 0. The remaining denominator on line 5 in D3 is brought down to 1 
by {516}; j2 > 1 by (1 - j2){236} + j33+2-{326} + j66+2-{623}; and jg > 1 
by {156} - {534} + (2 + l-/m2){dim}. However, that may revive massless 
numerators ji^s < 0. To avoid infinite loops on recursive application of the 
relations, we eliminate single numerators by a general projection-operator 
method 

/dV d"p, /,|p;, (p, - qf] h[pl (P2 - 9)'l A^'iPuP^.q) = 

m^'Xz'i nl'"^' mM-.f] A^M, (3) 



where A{pi,p2i l) = ^Vi {dfiv — lulv I 'f) "92 1 ^'^d /i^2 are arbitrary functions 
of their scalar arguments. A numerator [{pi — P2YY t)e re-expanded: 

{pi - P2f = \{pI + {Pi - q? +pI + {P2 - qf 

- [pI - (pi - q)V2 - iP2 - qf]/q^ -q^- A{pi,p2,q)}. (4) 

Odd powers of A{pi,p2,q) fall out after integration, and for even powers 
Eq. @ yields A{p,p, q) = 2[p^ + [p — qY] — [p^ — {p — qYf / q^ — q^ ■ For we 
identify the left-hand side of Eq. (^) with ci, and q^ with C3. If we deal with 
ji = —1, js = 0, = Jq, then on the right-hand side of Eq. (|D generates 
the original integral which can then be eliminated. 

The denominator > 1 is reduced by {156} — {236} while ji < 0. At 
ji = with j2,4-6 = 1 we apply {236}, eliminate ji = —1 by Eqs. and (Q), 
reduce je = 2 as usual, and expand the first numerator again. The original 
integral with the same value of js can be expressed by solving the resulting 
equation. 

The case .03(0, 1, 1, 1, 1, 1) can be transformed as follows. The differential 
equation for the two-loop subgraph without line 2, Eq. (43) of Ref. H], is 
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divided by q^{q^ + im?) and integrated over q. Derivatives with respect to m? 
are taken via {dim}, and after substituting simple integrals we get 
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D3(0, 1, 1, 1, 1, 1) = %^ D,{0, 1, 0, 1, 1, 1) 

2m 

/ ,x 3Af/2-5 

8[m'j r(|A^ - 1) r(5 - A^) 



(A^ - 2)(A^ - 3)(A^ - Af L r(3 - ^A^) 



(5) 



where each loop integral was divided by 7r^/^r(3 — ^N); in one loop, this 
modification agrees with the standard MS definition but is more conve- 
nient in higher-loop massive calculations. As a result of the transformations, 
any Da-type diagram is reduced to simpler types and two master integrals 
^3(0, 1,0, 1,1,1) and 1^3(1,. ..,1). 

With j3_4 < 0, i?4 is a particular case oi Bm- The numerator on line 5 can 
be taken off by {120} as j2 > 1 (massless lineO with index is assumed to 
connect 1.3 and 2.4); by {210} as ji 7^ 1; by {430} as ^3 > 1; or {dim} should 
be applied otherwise. The ji < numerator is eventually reduced to ji = 
(hence, to the two-loop massive bubbles) by {120} if j2 > 1, or by {210}. The 
denominators ji,2 > 1 can be brought down to 1 by solving {120} relative to 
2~^5~ , or {210} — 2m^5^{120} relative to 1"'"5~. That always increases j'5. 
The latter helps to reduce the denominators in the massive one-loop subgraph 
by 2{340} — {430}. The extra denominator > can be integrated off by 
applying {dim} and reducing j2 = 2 by {120}. Iteratively, we arrive at the 
master integral £"4(1, 1, 1, 1, 0) = 0^(0, 1, 0, 1, 1, 1). 

Further prototypes are quite simple indeed. The numerators are manage- 
able by adjacent-triangle relations, and the denominators can be reduced by 
combining three relations for a triangle that contains the line. The master 
integrals are (1, •••5 1), E^^l, and Dj^{l, In D2 the massless 

relation {125} solved with respect to the free term allows one to cancel out 
a denominator. Thus, in the end D2 is reduced to F functions, just as Di 
does. For E2 a massless relation can be constructed as {524} — {dim}. 

The described algorithms are implemented as a package of procedures in 
the symbolic-manipulation language FORM well suited to evaluating the 
Feynman diagrams as well as any polynomial-like expressions with a large 
number of terms. However, the efficiency of the essentially recursive programs 
is restricted by some features of the existing FORM translator. In particular, 
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'infinitely' iterative substitutions are only allowed without any intermediate 
sorting of terms. On the other hand, the recurrence relations generate rather 
many equal terms, and after exceeding certain machine- dependent limits on 
the size of the scratch expression generated in a module, the sorting becomes 
extremely slow. The only way out is to use step-by-step sorting inside the 
preprocessor T^do loops with a pre-estimated number of repetitions. But 
sometimes the number is rather difficult to guess at beforehand, while any 
misjudgement spoils the program performance. 

Therefore, it would be highly desirable to implement a kind of the prepro- 
cessor 9^repeat/endrepeat construct into FORM, which would terminate 
as soon as no actual transformations are applicable in any module inside its 
bodyQ Also of use would be any means of redefining preprocessor variables, 
based on global tests on the terms of sorted expressions. An invariable es- 
sential inconvenience for structured packages is the global scope of all names 
in FORM. 

The first application of the described package was the evaluation of the 
three-loop QCD correction to the electroweak p parameter |^: 

^QCD = + 2 C(2)] ^+{m-m C(2) - W C(3) + W C(4) 
- |C(2) ln2 + Ml 5^ - |i?4 - 3|/^3 - [i^ - f C(2) + |C(3)K 

- - I ^/) [1 + 2 C(2)] In (/xVm?) } ( ^ ) (6) 

where a<j is the QCD coupling constant at the renormalization scale /i in the 
MS scheme with the total number of quark fiavors nj (=6); is the pole 
mass of the top quark; B^^ has been introduced in Ref. 0; ^2 determines 
the finite part of the two-loop massive bubble master integral ; and is 
the finite part of 1)3(1, 1, 1, 1, 1, 1) which has been evaluated numerically by 
the momentum-expansion method |I[ and independently in Ref. |p. An error 
in the coefficient of C(4) in the original publication has been fixed, so that 

I thank Timo van Ritbergen for informing mc that in an wxperimental version of 
FORM 2.2 it is possible to terminate #do loops as nothing changes. However, this 
undocumented feature is unavailavle in public versions and, as Jos A.M. Vermaseren com- 
municates, liable to changes in FORM 3 
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Eq. completely agrees with the independent calculation of Ref . |Q . 
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